Master equations for pulsed magnetic fields: Application to magnetic molecules 
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We extend spin-lattice relaxation theory to incorporate the use of pulsed magnetic fields for 
probing the hysteresis effects and magnetization steps and plateaus exhibited, at low temperatures, 
by the dynamical magnetization of magnetic molecules. The main assumption made is that the 
lattice degrees of freedom equilibrate in times much shorter than both the experimental time scale 
(determined by the sweep rate) and the typical spin-lattice relaxation time. We first consider the 
isotropic case (a magnetic molecule with a ground state of spin S well separated from the excited 
levels and also the general isotropic Heisenberg Hamiltonian where all energy levels are relevant) and 
then we include small off-diagonal terms in the spin Hamiltonian to take into account the Landau- 
Zener-Stiickelberg (LZS) effect. In the first case, and for an S = 1/2 magnetic molecule we arrive at 
the generalized Bloch equation recently used for the magnetic molecule {Ve} in Phys. Rev. Lett. 
94, 147204 (2005). An analogous equation is derived for the magnetization, at low temperatures, 
of antiferromagnetic ring systems. The LZS effect is discussed for magnetic molecules with a low 
spin ground state, for which we arrive at a very convenient set of equations that take into account 
the combined effects of LZS and thermal transitions. In particular, these equations explain the 
deviation from exact magnetization reversal at B w observed in {Ve}. They also account for the 
small magnetization plateaus ("magnetic Foehn effect"), following the LZS steps, that have been 
observed in several magnetic molecules. Finally, we discuss the role of the Phonon Bottleneck effect 
at low temperatures and specifically we indicate how this can give rise to a pronounced Foehn effect. 

PACS numbers: 75.45.+j, 75.50.Xx, 75.60.Ej 
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I. INTRODUCTION 



The subject of magnetic molecules has attracted much 
attention both for its scientific importance for studying 
fundamental issues in nanomagnetism as well as for po- 
tential applications. Within each molecular unit are em- 
bedded a finite number of magnetic ions, coupled via 
Heisenberg super-exchange interactions. Furthermore, 
the intermolecular magnetic interactions are of dipolar 
origin and can usually be neglected. As a result, measure- 
ments on crystalline samples reflect the magnetic proper- 
ties of isolated individual molecules, with its most promi- 
nent feature, arising from the finite number of magnetic 
ions, being the appearance of a discrete magnetic energy 
level spectrum. 

This feature of the spectrum is reflected in the relax- 
ational behavior which mainly arises from the interac- 
tion with environmental degrees of freedom, i.e., a "heat 
bath", such as phonons: The relaxation times of the 
dynamical magnetization can become very long even at 
moderately high temperatures. Specifically, in molecules 
with a high spin ground stat o 1 ^ 3 (single molecule mag- 
nets (SMM's)) an anisotropy energy barrier is responsi- 
ble for relaxation times as high as 10 3 — 10 5 sec, whereas 
in some molecules with a low spin ground state, they 
are of the order of 1CT 3 - 10° sec for T < 4 K. 4 i 5 i 6 i 7 i 8 i 9 
The existence of long relaxation times becomes mani- 
fest through the appearance of dynamical hysteresis ef- 
fects when using pulsed magnetic fields, and this is one 
of the most exciting phenomena observed in magnetic 
molecules. Clearly, the hysteretic behavior is observable 
when the experimental time scale r e (determined by the 



field sweep rate) is in the regime of the spin-lattice re- 
laxation times (or shorter) . This opportunity is available 
for magnetic molecules in conjunction with the current 
experimental capability of using strong magnetic fields 
with sweep rates as high as 1 Tesla/ms^iSiS 

Another effect manifested in pulsed field measure- 
ments, is the appearance of abrupt magnetization steps 
at given fields pLSiiLLS^ These steps are quantum- 
mechanical in origin and reveal the existence of small, 
off-diagonal terms in the spin Hamiltonian which in 
turn give rise to avoided level crossings in the mag- 
netic energy spectrum and to Landau-Zcner-Stuckelberg 
(LZS^&ii transitions. The origin of these interactions 
may be single-ion anisotropy or anisotropic exchange. 
The characteristic energy splitting 5 of the LZS effect 
varies greatly among magnetic molecules. For istance, in 
SMM's due to the large spin of the ground state, S/ks 
can be of order 1CP 7 K (as usual, fee stands for Boltz- 
mann's constant), whereas in molecules with low spin 
ground state, typically S/ks ~ 0.1 Km& This implies 
that for currently available sweep rates one can probe the 
non-adiabatic regime in SMM'sji^ whereas in molecules 
with low spin ground states we are already in the extreme 
adiabatic regime. Experimentally, and for the molecules 
with low spin ground states, there have been observed 
deviations from the pure quantum-mechanical prediction 
regarding the height of the steps that have been associ- 
ated with the role of dissipation whithin the LZS regime. 

Of particular interest is also another effect induced 
by dynamic fields, namely the appearance of small 
plateaus^ai in M(t) following each magnetization step, 
thus giving rise to satellite peaks in dM/dB. It has 
been first showni for the low temperature experiments 
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on the magnetic molecule { V15} that the origin of these 
plateaus is the Phonon Bottleneck (PB) effecti 4 i 5 i 13 i 14 
Numerically solving a quantum master equation that 
had been previously derived for static fields, Saito and 
Miyashita 15 provided an alternative viewpoint of this ef- 
fect, which they termed "the magnetic Foehn effect": 
They suggested that this behavior is widespread whether 
or not one is in the PB regime and that it is a conse- 
quence (or an "after-effect") of the LZS transitions. As 
it turns out (see Sec. IV below), the PB effect can give 
rise to an enhancement of the Foehn effect. This has 
been recently observed experimentally and will be re- 
ported elsewhere^ 

At present, a first-principles account of such relax- 
ational phenomena in magnetic molecules induced by dy- 
namical magnetic fields is lacking. Our main goal is to 
show that one can generalize the conventional spin-lattice 
relaxation theory in the context of pulsed fields studies of 
magnetic molecules. The present work is devoted to (and 
motivated by) pulsed field studies of molecules with a low 
spin ground state: The simplicity of these systems, apart 
from providing a basis for better understanding the main 
physical ideas, allows one to directly compare the predic- 
tions of the generalized theory with experimental data. 
Hence, these systems, when subject to pulsed fields, pro- 
vide a convenient means for obtaining information on the 
various relaxational processes and microscopic interac- 
tions present in all nanomagnetic systems. 

One such system is the (S = 1/2) magnetic molecule 
{V6}r*i£ which shows both pronounced hysteresis loops 
as well as nearly complete reversals of the magnetiza- 
tion at B w 0. The hysteresis loops have been accurately 
reproduced 6 using a generalization of the standard 1318-19 
Bloch equation which in turn revealed that the one- 
phonon acoustic process is the dominant relaxation mech- 
anism at low temperatures, and in addition provided an 
estimate of the spin-phonon coupling energy. The first- 
principles derivation of this equation is provided here 
within the more general context of our analysis. On the 
other hand, the abrupt magnetization reversals at B w 
were interpreted as the result of adiabatic LZS transi- 
tions originating from the existence of a small (~ 0.4 K) 
intra-molecular, anisotropic exchange. The small devi- 
ation from the pure quantum-mechanical prediction of 
complete magnetization reversal was attributed to dissi- 
pation effects inside the LZS regime, but no quantitative 
account was given in Ref. |6:. This effect is also analysed 
in the present work. 

Another class of molecules with low spin ground state 
where the generalized theory can be easily applied is that 
of the antiferromagnetic (AFM) ring systems at low T. 
These are magnetic molecules comprising an even num- 
ber of uniformly spaced met al ions arrange d as a planar 
ring (see for example Refs. I20l2ll22l23l2^ . The AFM 
exchange interactions give rise, to a non-magnetic 5 = 
ground state, a first excited S = 1 triplet state, etc. 
In addition to their hysteretic behavior, these systems 
can show several magnetization steps and sometimes the 



small magnetization plateaus mentioned &boveM& As we 
show below, the present work accounts for these dynam- 
ical effects in a general way. 

The organization of this paper is the following. In 
Sec. II we develop the spin-lattice relaxation theory in 
pulsed fields for magnetic molecules with a spin S ground 
state that is well separated from the excited levels. We 
arrive at a generalization of the standard master equa- 
tions and show how these lead to the generalized Bloch 
equation for the case of S = 1/2 mentioned above. In 
Sec. Ill we extend this theory to include the general 
isotropic Heisenberg Hamiltonian where all energy lev- 
els are relevant. We apply this to AFM rings at low T 
and for fields around the first level-crossing field value, 
where they behave as two-level systems. We provide the 
treatment of dissipative LZS transitions in Sec. IV. This 
is done for the case of a level anti-crossing between two 
levels with different magnetic quantum numbers. We ap- 
ply the resulting theory to the spin 1/2 case and that of 
AFM rings at low T, and demonstrate how this theory 
accounts for the deviation from the quantum-mechanical 
prediction for the magnetization steps as well as the for- 
mation of the plateaus mentioned above (Foehn effect). 
In this way, in particular, we provide an explanation for 
the deviation from exact magnetization reversal observed 
in {V e }&We also indicate the role of the PB effect at low 
T, and specifically how it can give rise to a pronounced 
Foehn effect. Finally, in Sec. V, we provide a brief dis- 
cussion and summary of the present work. 



II. MASTER EQUATIONS FOR SPIN S 

We first consider a magnetic molecule with a ground 
state of definite total spin S, well separated from the 
excited levels^ We assume that the temperature is low 
enough and the field regime covered is well below the low- 
est level-crossing field, so that we only need to consider 
the ground state spin S level. We focus on the isotropic 
case i.e., we do not consider non-diagonal terms in the 
spin Hamiltonian. 

We employ the standardSSiSSiSl method of treating 
both the spin and the bath degrees of freedom quantum- 
mechanically. The Hamiltonian of the combined system 
(spin S + heat bath) is written as 

H(t) = H s (t) + H B +V , (1) 

where H s (t) = h^B(t)S z = hf(t)S z corresponds to the 
Zeeman energy, Hb is the bath Hamiltonian, and V 
is the spin-bath coupling. As usual, 7 = gfis/^ de- 
notes the electronic gyromagnetic ratio. In our notation, 
f(t) = 7-B(f) has units of frequency and the spin opera- 
tors are taken to be dimcnsionlcss. Typical sweep forms 
are shown in figures below; the experimental sweep time 
r e can be as short as 1 ms. 

For our purposes it is unnecessary to specify the de- 
tailed form of the interaction V, however, it can always 
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be written in the general form 

v = hJ2 A q ®R<i, ( 2 ) 

where A q (R q ) are hermitian operators of the spin (bath) 
system. This coupling may, for example, originate from 
the modulation of the exchange coupling terms (between 
the individual magnetic ions of the molecule) or the mod- 
ulation of the interaction of each individual ion with its 
local environment (e.g. the crystal field) which dynam- 
ically affects the spin degrees of freedom through the 
spin-orbit coupling. 14,28 In both cases, the modulation 
originates from lattice deformations, i.e., phonons. The 
terms bath and environmental degrees of freedom will be 
understood to mean phonons. 

The equation of motion for the density matrix ptot(t) 
of the combined system follows the von Neumann 
equation 2 ^ 2 ! 

Ptot{t) = -i[H{t)/n,ptot{t)] ■ (3) 
We switch to the interaction picture 

p tot = e Hnt)S.+H Bt /K) ptQt e -«F(t)S.+H B t/h) t (4) 

where F(t) = J* f(t') dt' . The equation satisfied by ptot 
is 

P tot (t) = -i[V(t)/h,p tot (t)], (5) 
where we have defined 

V{t) = hJ2A q (t)®R q (t), (6) 

with A q (t) = e lF ^ s * A q e ~ lF ^ s % and similarly R q (t) = 
e iH B t/h e -iH B t/h^ n 0Wj a u tn e relevant information 
about the spin system is contained in the so-called re- 
duced density matrix p(t) = Tr;, (ptot{t)), since for ex- 
ample 

< S z >= Ti- s Tr b {p tot (t)S z } = Tv s {p(t)S z } , (7) 

where Tr s (Trf,), denotes the partial trace over the spin 
(bath) degrees of freedom. Thus we are mainly interested 
in finding an equation of motion for pit). We first make 
the assumption that the spin-bath coupling is sufficiently 
weak so that the Born approximation can be used. Fur- 
thermore, we consider temperatures high enough (T > 1 
K) so that the number of available phonon modes per 
spin is very large and as a result the phonons equilibrate 
independently from the spins. Given this consideration 
and the expectation that the phonon relaxation times T& 
(typically Tf, < 10 -6 sec) are much shorter than both the 
experimental time scale r e and the spin-lattice relaxation 
time t s , the density matrix of the compined system (spin 
system + bath) can be factored as ptot{t) ~ p(t) <8> Pb, 
where ps — e~@ HB / Zb describes the stationary state of 
the heat bath at temperature T. Here Zb denotes the 



bath partition function, and (3 = l/(fcsT). The above 
factorization of the total density matrix is expected to 
break down at sufficiently low temperatures (typically 
T < 1 K) where one expects the phonon bottleneck (PB) 
effect 413-14 to take place. Employing the above approx- 
imations, one arrives at the following integro-differential 
equation of motion for p^^SL 

1 f* _ 
p(f)= duTr b [[V(u),p(u)®p B },V(t)} . (8) 

h Jo 

We denote the adiabatic eigenvalues of H s (t) by e m (t) = 
hf(t)M, where M = —S, . . . , S. The adiabatic excita- 
tion frequencies of H s (t) are of the form oj^t) = f(t)p, 
where p = —2S, . . . , 25. Now, for a given operator A a 
of the spin system it is very convenient to construct the 
so-called eigenoperator A a ^ corresponding to a given ex- 
citation frequency w^(t) as 

A q ^ = {A q ) M M>\M >< M'\ 6m>-m,,> , (9) 

M.M' 

where Sjj denotes the Kronecker delta symbol. These 
operators obey the equation 

^, = E^, (io) 

n 

where the sum extends over all possible integers p. The 
reason for introducing the eigenoperators A q ^ is that 
they take a very simple form in our interaction picture, 
namely 

A q ^(t) = e- iP ^A q!fl . (11) 

We note in particular that for a static external field £>o, 
we have F(t) — uj^t, where = jBq and the phase 
factor in Eq. I|ll|) becomes the familiar form exp(— iu>ot). 
In the present case, Eq. JSJ along with Eqs. lfTU|) 
and Ijllfl and the variable change u — > t — u, give 

p{t) = V / due-^W"' e -«"(*-« (A q , ifi p(t- u)A q ^ 
i i Jo 

-A q ^,A q ,^p{t - u)) < R q (u)R q ,{0) > B +/i.c(I2) 

where the quantities in angular brackets 

< R q (u)R q ,{0) > s = Tr b (p B R q (u)R q ,{0)) , (13) 

are equilibrium time correlation functions of the bath, 
and the symbol h.c. denotes hermitian conjugate. As 
mentioned already, Tf, << t s ,t c . This allows one to per- 
form the following simplifications. First, we may extend 
the upper limit of integration in Eq. (|12f> to infinity. 
Second, the variation of p(t — u) in the time scale of Tf, 
is extremely small (since r s >> u), so it is justified to 
perform the usual Markov approximationjS^SiSi namely 
we replace p(t — u) by p{t). Similarly, the variation of 
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F(t — u) in the time scale of t& is also small. Therefore 
we may approximate F(t — u) by 

F(t - u) « F(t) - u ■ F(t) = F(t) - uf(t) , (14) 

i.e., by a Taylor expansion through first order in u. The 
neglect of higher order terms is valid since u 2 f(t) « 
uf(t), is equivalent to Tb << r e (with Tb taken as the 
maximum value of u, and r e as a typical value of ///). 
Substituting Eq. (jTlJl into Eq. (jT2|) . we obtain 

j>(t) = J2 e- iF «^r gg ,(^))x 

qq',HH' 

(A q , 4l p(t)A q . y - A q ^,A g , tl ,p(t)) + h.c. , (15) 

where we define the bath correlation functions, at the 
time-dependent frequency o;^(t), as 

/>oo 

1X^(0) = / du u < i? 9 (u)i? 9 ,(0) >s • 

Jo 

(16) 

We emphasize that the 'adiabatic' uj^t) factor in Eq. 
I|lfi|l originates from the second term in the Taylor ex- 
pansion of F(t — u) and it has important implications in 
what follows. In particular, since is proportional to 

the instantaneous field B(t), it will lead to an equation 
of motion with relaxation rates that depend explicitly on 
Bit). 

We further adopt the so-called rotating wave approx- 
imation (RWA)i 25 i 26 i 27 The relative rate of change of a 
typical phase factor of Eq. (fl5|l . with p, + p! ^ is pro- 
portional to 7-B(i) which is of order 10 11 s _1 (for B ~ 1 
Tesla). This implies that such non-secular terms "os- 
cillate" very rapidly during the experimental time scale, 
and thus do not appreciably contribute to the dynamics 
in Eq. ((15(1 . 33 Thus, by retaining only the terms with 
p! = —p, we may replace Eq. 1(15(1 by 

11' ,M 

-A q ,- M A g , t pp(t)) + h.c . (17) 

This is our generalized master equation in the weak cou- 
pling and RWA limit and for slowly changing (compared 
to Tb) external fields. 

Before we discuss the major consequences of this gener- 
alized master equation, we go one step further and derive 
the equations of motion for the populations p mm of the 
various states \M >. Using the matrix elements 

{Aq,n)MM> = {A q )MM'^M'-M,ij,, (18) 

and the relation pmm = Pmm, one finds that the popu- 
lations decouple from the non-diagonal terms and evolve 
according to the following generalization of the standard 
Pauli master equation 

Pmm = W M 'M(t)pM'M> - Wmm' if) Pmm ■ (19) 

A-/' M' 



The transition rates Wm->M' = Wmm> are defined as 

WmM> = 7gg' (UMM> (t)){A q ) M M> (A q >) M 'M, (20) 

1,1' 

where Uwmm' = e M — cap, and 

7gg'M = T qql (uj) + r* lq {uj) = 

/ + 0O 
due luu < R q (u)R q ,(0) > B . (21) 
-oo 

An important property of the transition rates is the de- 
tailed balance condition that arises from the following 
quantum property^SSiSl of the (stationary) bath corre- 
lation functions appearing in Eq. I(21|l 

Jqq'(-oj) = exp(-/3?iw)7 < j/ g (w) . (22) 

The detailed balance condition follows straightforwardly 
from Eq. J2jOJ alon S with Ec l- <G2> 

W M 'M = exp(-f3TwjMM>)W M M> ■ (23) 

We have now arrived at the main results of this Section. 
A comparison with the standard relaxation theory (for 
static fields) shows clearly the physics underlying the gen- 
eralization made here: Equation (121 )[) involves transition 
rates that depend on the adiabatic energy excitations of 
the spin system, which in turn are proportional to B(t). 
All information about specific details of relaxation mech- 
anisms is contained in the bath correlation functions as 
well as the matrix elements (A q )MM' (see Eq. I|2()|l ). 

Reviewing our derivation, the key assumption made 
is that the bath degrees of freedom equilibrate indepen- 
dently from the spins in times 77, << r e ,T s . The basis 
of this assumption is that the number of available envi- 
ronmental modes per spin is so large that the bath can 
be considered to be a large reservoir, hence allowing the 
neglect of any feedback from the spin dynamics. This 
allowed the decomposition ptot(t) « p(t) (g> pb and also 
the Taylor expansion (particularly the second term) in 
Eq. 1(14(1 . Despite the intuitive appeal of this assump- 
tion, it must be checked by comparing the predictions of 
the resulting theory with experimental data. 

A simple, realistic system where the present theory is 
easily applicable and is in fact in excellent agreement 
with experimental data (for T > 1.5 K) is the magnetic 
molecule {V 6 }£ mentioned in the Introduction. It is 
straightforward to show that the generalized Bloch equa- 
tion used in Ref. 0, to reproduce the experimental data 
for this S = 1/2 system follows immediately from the 
present theoretical framework: For S — 1/2, the equation 
of motion (Eq. iJUJl) for the populations of the spin- up 
(|+ >) and spin-down (|— >) states reads 

P++ = W-+(t)p— - W + -(t)p++ , 

P— = W+-(t)p ++ - W-+(t) P — . (24) 

On using the normalization condition p ++ +p = 1, and 

the detailed balance condition W-\ = exp(/3ftfi)W |_ 
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(here = 2^B) one finds that the magnetic moment per 
spin, M = —Tij < S z >, follows the equation 

M(t) = Ts{T l iB{t)) ( M e q (T, B{t)) - M{t)) , (25) 

where M eq = (hj/2) tanh (/3M2/2) and and l/r s = 
W- + (I + cxp(f3HQ)). This is the generalized Bloch 
equation that was used in Ref. 0. Its physical interpreta- 
tion is that M(t) relaxes towards the instantaneous equi- 
librium value M eq (T, B(t)) with a relaxation rate that 
depends explicitly on B(t). 

In principle one can invert Eq. (|25|l and extract 1 /t s in 
terms of M{t) and M(t) obtained by experiment and the 
adiabatic equilibrium magnetization M eq (T, B(t)). Al- 
ternatively, one can directly compare the experimental 
data with a numerical solution of Eq. (|25|l by choosing 
a physically appropriate functional form of l/r s (T, B{t)) 
and adjusting the free parameters. Due to the explicit 
dependence of l/r s and M eq on B(t) one can obtain in- 
formation on the underlying specific relaxation mecha- 
nism^) by using different sweep forms B(t). Along these 
lines, it was confirmed in Ref. |6j that for the magnetic 
molecule {Ve} and for 1 < T < 5 K, the dominant con- 
tribution to 1/t s is the one-phonon processes term 1/r* 
given by 

1/r* = AQ{t) 3 coth((3hn(t)/2) , (26) 

with A = 3Vl l /(2nhpv 5 ), where v denotes the sound 
velocity, p the mass density, and V s ; the characteris- 
tic energy modulation of the given spin-phonon cou- 
pling mechanism™*^ Apart from the establishment of 
the dominant relaxation mechanism at 1 < T < 5 K, 
a first estimate of V s ; (~ 0.35 K) was obtained. More 
generally, the excellent agreement of this theory with ex- 
perimental data signifies that our starting assumptions 
are valid for {Ve} for T > IK. However, our main as- 
sumption, namely that the phonons remain in equilib- 
rium at all experimental times, can be expected to break 
down at lower T since the number of available resonant 
phonons per molecule rapidly decreases on cooling, and 
the phonon bottleneck (PB) effedpiiiii takes place. In 
fact, preliminary data™ for {Ve} at T = 0.6 K shows a 
significant deviation from the theory suggesting the onset 
of the PB effect. 

We remark that the theory of this Section, which is 
based on the Hamiltonian of Eq. Q), cannot account for 
the magnetization steps observed in {Ve} at the level- 
crossing regime B w 0. These, however, can in fact be 
explained in terms of adiabatic LZS transitions. The nec- 
essary extension of our theory is given in Sec. IV. 

III. MASTER EQUATIONS FOR THE 
GENERAL ISOTROPIC HEISENBERG MODEL 

Here we extend the previous analysis and discuss the 
general isotropic Heisenberg model where all the energy 



levels are relevant. The analysis is parallel to the above 
and straightforward, and thus only the main new points 
are emphasized. 

The Hamiltonian of the combined system (magnetic 
molecule + heat bath) is again given by Eq. Q where 
now H s (t) explicitly includes the Heisenberg exchange 
Hamiltonian 

H = JijSi ■ Sj , (27) 
i<j 

and where Jy denotes the exchange constants between 
the spins at sites i and j. The eigenstates \n > of H s 
are of the form \n >= \v S M >, where v corresponds to 
additional quantum numbers. Since Hq and S z commute 
one can define the interaction picture by 

p- e i(H t/H+F(t)S,) p e -i(H t/h+F(t)Sz) ^ ^gs 

In addition, the adiabatic energy eigenvalues have the 
form e n (t) = eon + hf(t)M, where the first term corre- 
sponds to the zero-field spectrum of the exchange Hamil- 
tonian (Eq. HZl) and the second term to the Zee- 
man splitting energy. Hence, the adiabatic excitation 
frequencies are of the form uj(t) = luq + f{t)ji, where 
Jiluq = eon' ~ e 0n and A 1 = M' — M . Thus, any given exci- 
tation frequency Lu(t) can be characterized completely by 
<x>o and p. As before, we introduce a set of eigenoperator 
A q (uo,fi<) given by 

A q (uj ,p) = y^ j { A i)nn'\n >< n'\ x 

n,n' 

hl'-M.fj. 8( e 0n> -£On,ftwo) , (29) 

which take the following form in the interaction picture, 
A q (u Q , p) = e- iuot e- iF ^A q (u;o, p). (30) 

Employing the same steps as in Sec. II, one arrives at 
the master equation 

<?<?' oj{t) 

-A q (-Qj , -p)A q ,(uj ,p)p(t)) + h.c.(31) 
Then, using the matrix elements 

(A q (LO , p)) nn i = [A q ) nn ,8M'-M,n^{ e W ~ <^On,flUJo) (32) 

one obtains the generalized Pauli equations as before 

Pnn = J]] W n 'n(t)Pn'n' ~ W nn'(t)Pnn , (33) 
n' n' 

where W nn ' are given by 

W nn > = y^7gg' (^W {t))(A q ) nn i (A?')n'n , (34) 
qq ' 

with 7 gg /(w) as in Eq. I|21|) . According to the above 
analysis, the dynamics of the reduced density matrix has 
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FIG. 1: (Color online) Hysteretic behavior of AFM ring sys- 
tems at low T and for two commonly used sweep forms B(t) 
shown in Fig. 1(b). The loops are obtained by a numerical 
solution of Eq. 12511 . with M eq given in the text and 1/t s 
taken as the sum of the one-phonon term plus a small resid- 
ual contribution. All the parameters used are arbritrary. The 
equilibrium magnetization is denoted by the dashed line. The 
low-energy diagram and the true level crossing between the 
two relevant states |0, > and |1, —1 > is shown in Fig. 1(c). 



the same major features as in the case of Section II, but 
now with the appropriate and reasonable modification for 
the excitation frequencies. 

As discussed in the Introduction, a class of magnetic 
molecules where the above analysis can be easily applied 
is that of AFM rings at low T and for fields in the vicin- 
ity of a given level crossing value. Specifically, we will 
assume T << Ao/^b, where Ao denotes the first, zero- 
field excitation energy and consider fields in the vicinity 
of the first level-crossing field (B c — Aq /h-f) only, where 
the singlet |0, > intersects the M = — 1 (or |1, — 1 >) 
level of the S = 1 triplet state (see Fig. 1(c)); a simi- 
lar analysis can be employed for fields in the vicinity of 
higher level crossings. For these temperatures and fields, 
the AFM rings behave as a two-level system. Hence, 
we are dealing with a situation that is very similar to 
the spin 1/2 case, discussed in Sec. II. In fact, one can 
arrive at the same generalized Bloch equation for M(t) 
(Eq. l|55)l). where now M eq = Try sech(/3A(i)/2), and 
A(t) = Ao — hjB(t). Similarly, l/r s will depend explic- 
itly on A(t). For istance, the contribution 1/r* of the 
one-phonon processes to the relaxation rate is given by 
Eq. (H)J) with n(t) replaced by A(t)/h. 

Similarly with the approach followed in Ref. |(| for 
{V6J-, finding the physically appropriate functional form 
for 1/t s can be facilitated by comparing the theoretical 
predictions with experimental data for a variety of field 
sweeps. We illustrate this idea in Fig. 1(a) by show- 



ing typical hysteresis loops obtained by numerically solv- 
ing Eq. i|25|) using two different, commonly used sweeps 
shown in Fig. 1(b). The relaxation is assumed to be 
driven by one-phonon processes plus a small residual 
term. Since we have assumed an isotropic Hamiltonian 
and therefore a true level-crossing at B rj B c , nothing 
exciting happens in that regime. Instead, our choice for 
1/t s gives rise to a plateau for B « B c since in that 
regime the one-phonon term vanishes and the relaxation 
is driven only by the small residual term. This is quali- 
tatively different from observations in AFM ring systems 
such as the dimer {Fe2}^ and dodecanuclear {Fei2}^ 
which show a step-like behavior of M(t) and consequently 
peaks in dM(t)/dB at the level crossing fields. In addi- 
tion, in some cases, the major peaks are accompanied 
with small satellite peaks; a small plateau in M(t) is 
formed after each step (Foehn effect). As we show in 
the following Section, both of these features, namely the 
magnetization steps and plateaus can be explained by 
extending our theory to include small off-diagonal terms 
in the spin Hamiltonian. 

IV. DISSIPATIVE LZS MODEL IN THE 
ADIABATIC REGIME 

We now include small off-diagonal terms in the spin 
Hamiltonian in order to take into account the combined 
effect of LZS and thermal transitions. We will consider 
the adiabatic regime: As mentioned in the Introduction, 
this is indeed the relevant regime for a large class of 
molecules with low spin ground state, where the energy 
splittings are typically of order S/kB ~ 0.1 K. In partic- 
ular, the magnetic molecules {Vg}^ and { Vis)^ as well 
as the AFM rings fall in this category. Since the effect 
of the small non-diagonal terms becomes manifest only 
in the immediate vicinity of the intersection of twoM en- 
ergy levels, we are able to construct the following quite 
adequate theory. 

A. General theory 

We consider two energy eigenstates of the isotropic 
Heisenberg Hamiltonian (denoted as \mi > and \m,2 >) 
with total spins Si and £2 and magnetic quantum num- 
bers mi and 777*2 , respectively, which are coupled by a 
small off-diagonal (anisotropic) term. This term gives 
rise to an avoided level crossing between the two energy 
levels with a small energy gap denoted by S. For fields in 
the immediate vicinity of this level anti-crossing, one can 
write the Hamiltonian in the basis of \m\ > and \rri2 >, 
i.e., 

where Eq denotes the zero-field energy difference between 
the two states for 5 = 0. Without loss of generality, 
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we assume that m% > mi, with Eq,S being real and 
positive. In the absence of the off-diagonal term (i.e., 
(5 = 0) the two levels cross at the moment when / = 
f c = Eo/h(m2 — mi). We denote the adiabatic energy 
levels of H s (t) by e± (t) and the corresponding eigenstates 
by \e±;t >, i.e., H s (t)\e±;t >— e±(t)\e±;t >. In the 
basis of |mi > and |m 2 >, they can be expressed in the 
convenient parametric form 

. . _ / cos 9/2 \ 
\ e +' t> - [ sin 9/2 ) ' 

i^>=("c»//2 2 )' (36) 

and e±(t) = [Eq + hf(mi + 7112) ± hd]/2, where 

n = [{5/fif + (mi - m 2 ) 2 (/ - / c ) 2 ] 1/2 . (37) 
The time-dependent parameter 9 is given by 

= ^1 S/ t f f w (38) 

(mi - m 2 )(/ - fc) 

and extends from to tt as / goes fron — 00 to +00. In 
particular, 9 w ir/2 when / is in the immediate vicinity 

Of fc. 

As in Sees. II and III, we switch to the interaction 
picture 

p = UjpU s , (39) 

where U s (t) is the evolution operator for the spin Hamil- 
tonian H s (t) alone, which obeys iKU s = H s (t)U s . Con- 
trary to the previous isotropic cases, it is clear that due 
to the presence of two non-commuting terms in H s (t) 
the form of U s (t) cannot be written in a closed analytical 
form.' 55 Thus, an analysis parallel to that of the previ- 
ous sections cannot be readily employed. Nevertheless, 
it is possible to circumvent this difficulty by exploiting 
the fact that we are in the adiabatic regime. According 
to the adiabatic theorem one has 

U a (t)\a;-oo >= e-^ (t) |a;< >, (40) 

where the phases 4><x(t) are given by <f> a (t) — 
J_ e a (t')dt' /%. As we mentioned before, adiabaticity 
holds even inside the LZS regime. 

We now express the spin operators A q appearing in 
the spin-phonon interaction term V (Eq. @), in the 
adiabatic basis as 

A i = J2 A f (t)\^t><0;t\, (41) 

where A^(t) =< a;t\A q \(3;t >. In the interaction pic- 
ture these take the form 

A q (t) = Af^e^-V-teW \a; -00 >< (3; -oo| . 

(42) 



Following the previous steps, it is straightforward to de- 
rive an expression similar to Eq. l|12|l , where one encoun- 
ters typical matrix elements such as A" & (t — u) and phase 
factors of the form e % ^ a ( t_tt ). The next step is to approx- 
imate A*P(t — u) » AgP(t), and also to make a Taylor 
expansion of the phases in first order in u, as before, i.e., 

<t) a {t-u)^(j3 a (t)-u€ a (t)/ri. (43) 

It is then possible to obtain the following set of 
Pauli master equations by introducing the representation 
p a p(t) =< a; — oo|/?(t)|/3; —00 >, and then performing 
the RWA, 

P aa = W Pa {t)pi3p - W af3 (t)p aa , (44) 
where the transition rates W a p are now given by 

W aP (t) =Y J Af(t)A f3 «(t) lqq ,{u Jaf} (t)), (45) 

qq> 

and jq q r(u)) as in Eq. 1|21|) . We now compare this ex- 
pression for the transition rates with the correspond- 
ing ones found previously. The new ingredient here is 
the extra time dependence of the transition rates car- 
ried by the matrix elements A q @(t). Clearly this results 
from the explicit time dependence of the adiabatic en- 
ergy states. Physically this means, for example, that 
longitudinal fluctuating fields (contained in V) can be- 
come effective in inducing transitions between the two 
levels inside the LZS regime because of the admixing of 
the two states. This introduces an additional compli- 
cation when one attempts to quantitatively account for 
magnetization data inside the LZS regime. 

It turns out that one can obtain a simplified Bloch type 

of equation for the quantity n(t) = p — /?++, in terms 

of which the magnetization can be simply expressed. Us- 
ing the above master equations and the normalization 

condition p + p ++ = 1, one arrives at the following 

equation for n(t) 

T s (t) 

where n eq = tanh(/3^/2), l/r s = W-+ (l + e (ihn ). In 
order to express < S z >= Tr{pS z } in terms of n(t) one 
notes that 

<S z >=J2<^t\p\(3;t> S^ a (t) 

a, (3 

= J2^ f,(t) ~ Mt)) p af 3S^ a {t) . (47) 

Now, for the longitudinal magnetization one can safely 
repeat the RWA by neglecting the terms with a ^ (3. 
Then < S z >~ p__S~~ + p ++ S+ + , or for the magneti- 
zation M, 

M a ^[(mi + m 2 ) - (m 2 - mifn(t) f --A] . (48) 
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plus a small residual term to account for the relaxation 
inside the LZS regime. The parameters chosen are some- 
what arbritrary; the exact shape of the M vs B curves is 
determined by the actual parameters of a given system 
(magnitude of 8, spin-phonon coupling terms, etc). As 
explained below, Eqs. I(46|l and (|48(l account nicely for 
all the dynamical effects shown in pulsed field measure- 
ments, namely hysteresis loops (already discussed in Sees. 
II and III) outside the LZS regime, the thermal deviation 
of the LZS steps from the pure quantum-mechanical pre- 
diction (see in particular Figs. 2 and 3), as well as the 
formation of magnetization plateaus (Foehn effect) im- 
mediately after exiting the LZS regime (see Figs. 2 and 
3, but mostly Fig. 4). Furthermore, we will discuss how 
the PB effect, which takes place at very low T, can give 
rise to an enhancement of the Foehn effect. 



FIG. 2: (Color online) Hysteresis and LZS effect at B ps 
for the case of an S=l/2 magnetic molecule, as obtained by 
solving numerically our equations (see text) for the full-cycle 
sweep shown in Fig. 2(b). The shaded area indicates the 
LZS regime. The equilibrium magnetization is denoted by the 
dashed line. Note the deviation from the exact magnetization 
reversal at B ~ as described in the text, and the formation 
of the small plateau after the step. B* denotes the field at 
which M crosses M eq , as discussed in the text. The low- 
energy diagram and the level anti-crossing between the two 
levels | - 1/2 > and | + 1/2 > is shown in Fig. 2(c). The 
small circles in Figs. 2(b) and 2(c) indicate the LZS regimes, 
where the magnetization steps take place. 



Equations 146JI and (|48|l are of central importance for 
describing the combined effects of LZS and thermal 
transitions. 36 According to Eq. I|48|l . the magnetiza- 
tion is given in terms of two distinct quantities, namely 
the ratio (/ — / c )/0 and n(t). Interestingly, the factor 
(/ — fc)/Q, carries the physics of the purely quantum 
LZS effect, since it changes sign inside the LZS regime as 
expected in the adiabatic regime of the LZS transitions. 
On the other hand, the quantity n(t) contains all infor- 
mation about thermal transitions and dissipation since 
its dynamics is determined by the spin-bath coupling ac- 
cording to Eq. firjfl. 



1. Hysteresis loops 

We begin by noting that for either 8 = or for fields 
outside the LZS regime one has (/ — / c )/f2 oc sgn(/ — / c ) 
and therefore the only time dependence of M stems from 
the quantity n(t). In addition, all matrix elements A^P 
appearing in Eq. (|45|l become time-independent. Then, 
by taking the time derivative of Eq. (|48|l . one recovers 
the results for the isotropic case, and in particular the 
generalized Bloch equation, Eq. I|25|l . derived before for 
the case of S = 1/2 and that of the AFM rings at low 
T. Thus, as expected, one can neglect the LZS effect for 
fields outside the immediate vicinity of level crossings. 
This justifies the use of Eq. (|25|l in Ref. 6, for fields away 
from B w 0. Typical hysteresis loops for fields outside the 
LZS regime are shown in Figs. 2(a) and 3(a) for the S = 
1/2 case and that of AFM rings at low T, respectively. 
One should note that although the magnetization obeys 
the same generalized Bloch equation (Eq. H25|l 1 f° r either 
8 = or for 8^0 but for fields outside the LZS regime, 
the solution is drastically different for these two cases 
(compare for example Figs. 1 and 3, for the case of AFM 
rings). This is because, when 8 ^ 0, the occurence of an 
LZS step introduces a different (as compared to the 8 = 
case) initial condition immediately after exiting the LZS 
regime. A direct consequence of this is the Foehn effect 
discussed below. 



B. Qualitative analysis 

We will now give a qualitative analysis of the predic- 
tions of this theory by examining the general structure of 
Eqs. Ij46|l and 1)48(1 . We will also demonstrate the main 
ideas by providing numerical solutions for typical exam- 
ples, namely the 5 = 1/2 case and that of AFM rings 
at low T. These are shown in Figs. 2, 4 and in Fig. 3, 
respectively. For these solutions, we have assumed that 
1/t s is the sum of the one-phonon process contribution 



2. Thermal corections to the LZS step 

A deviation from the exact quantum-mechanical pre- 
diction regarding the magnetization step at / ~ f c , is 
expected to arise from thermal transitions inside the LZS 
regime. This can be seen as follows. Assuming that one 
crosses f c from below, and denoting by Mi n and M out 
the magnetization when entering and when exiting the 
LZS regime, respectively, (similarly for rim and n out ), we 
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can obtain from Eq. I|48(l 

M in + Mout = -?vy[(mi + m 2 ) - (m 2 - mi) 2 5n LZ s /2] , 

(49) 

where Sulzs = n out — Tii„, denotes the overall change of 
n(t) inside the LZS regime, as obtained using Eq. (JJfjJ. 
This quantity is negative, i.e., n out < n in (this can be 
easily seen by plotting n eq vs / and solving Eq. 146(1 
graphically, i.e., without specifying the form of l/r s ). 
The first term of Eq. 1(491) gives the quantum-mechanical 
prediction in the adiabatic regime, since in the absence of 
thermal transitions inside the LZS regime (i.e., h(t) = 0) 
the second term vanishes. To be more specific, for the 
spin S = 1/2 case (ma = —mi = 1/2, and Eq = 0), Eq. 
gUJ gives 



1.0 



-M in - h~f\Sn L zs\/2 ■ 



(50) 



One then obtains the expected magnetization rever- 
sal {M ou t — —Mi n ) in the absence of thermal effects 
{Stilzs = 0). Thus, the second term of Eq. (|4*9"|) (or that 
of Eq. I(50(l) gives the thermal correction; its magnitude 
clearly depends on the competition between two time 
scales, namely r s and the time StLzs spent inside the LZS 
regime, given by St^zs = 2<y/(ft / yr), which is controlled 
by the sweep rate r. Hence, for a given T, the thermal 
correction becomes larger with decreasing sweep rates mm 
Furthermore, since l/r s is expected to increase with in- 
creasing T the thermal effects are more pronounced at 
higher T, as indeed observed for {V 6 }& Of course, at 
high enough T, the LZS effect is completely masked by 
the thermal transitions and the step disappears. More 
generally, it should be noted that it is Mi n + M out , rather 
than M ou t — Mi n (height of the step) , that is of more di- 
rect relevance in experimentally determining the extent 
of the thermal effects. The numerical solution for the 
5 = 1/2 case shown in Fig. 2(a) demonstrates the devi- 
ation from exact magnetization reversal at B rj 0, as a 
result of thermal transitions inside the LZS regime. This 
is consistent with the experimental data for {Vg}. 



3. The magnetic Foehn effect 

Another effect which is also of particular interest is the 
"magnetic Foehn effect"— mentioned in the Introduc- 
tion, which concerns the formation of plateaus shortly 
after exiting the LZS regime. Although these plateaus 
occur outside the LZS regime, they are a direct conce- 
quence of the LZS effect, as can easily be explained by 
considering the particular example shown in Fig. 4 (the 
analysis for the other cases shown in Figs. 2(a) and 3(a) 
is analogous). Immediately after exiting the LZS regime, 
M ut > M eq in contrast to the typical hysteretic behavior 
where M lags behind M eq . Then, outside the LZS regime 
where, as explained above, M tends towards M eq , M will 
decrease even though the field is swept towards larger val- 
ues. This drop of M eventually stops at a field B* , at 
the moment when M equals M eq . Thus, B* corresponds 
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FIG. 3: (Color online) Hysteresis and LZS effect at B w B c 
for AFM ring systems at low T and for fields around the 
first level-crossing value B c , as obtained by solving numeri- 
cally our equations (see text) for the half-cycle sweep shown 
in Fig. 3(b). The shaded area indicates the LZS regime. 
The equilibrium magnetization is denoted by the dotted line. 
Note the deviation from the pure quantum mechanical predic- 
tion regarding the step (see text) and the formation of small 
plateaus after each step. B* denotes the fields at which M 
crosses M eq , as discussed in the text. The low-energy dia- 
gram and the level anti-crossing between the levels |0,0 > 
and |1, —1 > is shown in Fig. 3(c). The small circles in Figs. 
3(b) and 3(c) indicate the LZS regimes, where the magneti- 
zation steps take place. This typical behavior of M(t) in the 
case of 8 7^ should be contrasted with that of S — shown 
in Fig. 1. 



to a local minimum in M (since then M = 0). Generally, 
this minimum will be observable if r s , in that regime, is 
much shorter than the sweep time r e . If this is not the 
case, a broad plateau around B* is formed instead. In 
fact, this seems (see below) to be the typical behavior 
in all experimental data reported so far. This behavior 
is shown in Fig. 4(a), where for a fixed T and with in- 
creasing sweep rate r, the minimum broadens with B* 
departing from B c = 0. One can exploit this feature ex- 
perimentally, in order to obtain an estimate for 5, since 
all B* corresponding to different sweep rates must lie on 
the equilibrium curve M eq (B(tj). 

The curves shown in Fig. 4 also suggest why a very 
steep minimum at B* (pronounced Foehn effect) has not 
been reported so far in realistic situations. In principle, 
as mentioned above, the minimum at B* should become 
more steep at lower sweep rates. On the other hand, 
on decreasing the sweep rate, one indirectly reduces the 
value of M ou tr~ thus being closer to M eq when exiting 
the LZS regime, and therefore compensating the effect of 
slowing down the sweep rate. 

Ideally, for the realization of a very steep minimum at 
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B* one needs (for a fixed sweep rate), a relaxation process 
that is slow enough at negative fields (i.e., large \M in \ and 
concequently large M out ) but becomes faster at positive 
fields, after exiting the LZS regime. Such an "asymme- 
try" in the field dependence of 1/t s cannot exist unless 
we abandon our main assumption that the phonons are 
in equilibrium at all experimental times. This is because, 
according to our analysis, with the phonons being in equi- 
librium, 1/t s depends on |-B(t)|, and thus must be sym- 
metric around B c — 0. This raises the question whether 
the PB effect, which takes place at very low T, could 
give rise to a pronounced Fochn effect. This is indeed 
supported by the analysis of the PB effect described in 
Ref. H for the magnetic molecule { V15}, which shows 
an asymmetry in the number of resonant phonons be- 
fore entering and after exiting the LZS regime (see, in 
particular, the inset of Fig. 2(b) of Ref. |j), as a re- 
sult of the strong coupling of the phonons to the spin 
degrees of freedom. In short, the physical origin of this 
asymmetry is that B c = is a level-crossing field for 
the energy spectrum of the spin system, whereas it cor- 
responds to a "reflection" point for that of the resonant 
phonons: Although, in general, when sweeping an ex- 
ternal field one progressively brings the spin system into 
resonance with different phonon modes (instantaneous 
resonance condition) this is not true in the vicinity of 
a level anti-crossing, since when exiting the LZS regime 
the relaxation is driven by the same phonons that were 
in thermal contact with the spins while entering the LZS 
regime (for more details, see Refs. 013). The above en- 
hancement of the Foehn effect due to the PB effect has 
been in fact recently observed experimentally and will be 
reported elsewhere^ 



V. SUMMARY 

We have extended the standard spin-lattice relaxation 
theory, in the context of pulsed field studies of magnetic 
molecules. Being easily applied to simple systems, this 
generalized theory can give important information on the 
underlying relaxation mechanisms and the microscopic 
interactions present in magnetic molecules in general. All 
the dynamical magnetization effects, including hysteresis 
loops, LZS transitions with or without dissipation and 
magnetization plateaus (Foehn effect), which are mani- 
fested in pulsed fields measurements, are accounted for 
by the comprehensive theory presented above. 

We first developed the theory for the isotropic case of 
molecules with a spin S ground state well separated from 
the excited levels but also the general Heisenberg model 
where all energy levels are relevant. We have shown, for 
two such simple cases, namely the spin S = 1/2 case and 
that of AFM ring systems at low T, that the dynami- 
cal magnetization obeys a generalization of the standard 
Bloch equation. In particular, this equation has been re- 
cently used£ for the magnetic molecule {Vg}, and was 
found to provide results in excellent agreement with ex- 
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FIG. 4: (Color online) Magnetic Foehn effect for a magnetic 
molecule with a spin 5 = 1/2 ground state, for fixed T and for 
three different fields sweeps with constant sweep rates r (in 
arbritrary units). The shaded area denotes the LZS regime. 
Because of the LZS effect, immediately after exiting the LZS 
regime (shaded area) M out > M eq , as opposed to a usual 
hysteresis situation. Then, outside the LZS regime, where 
M tends towards M eq (see text), M drops even though the 
field increases; a local minimum of M is formed at a field B* 
at which M = M eq . Note that this minimum is departing 
from the LZS regime of the magnetization step and broadens 
(eventually giving rise to a plateau) with increasing sweep 
rates. 

perimental data at T > 1 K, confirming that the domi- 
nant mechanism driving the relaxation is the one-phonon 
processes, and in addition providing an estimate of the 
spin-phonon coupling energy. Obtaining this information 
was greatly facilitated by using a variety of field sweep 
forms. 

We also extended the theory to include small off- 
diagonal terms in the spin Hamiltonian and thus take 
into account the combined effects of LZS and thermal 
transitions. This was done here for the large class of mag- 
netic molecules with a low spin ground state. For these 
molecules, and for the currently available sweep rates, 
one is in the extreme adiabatic regime. Our main inter- 
est in these systems, has been the role of dissipation on 
the LZS steps, as well as the formation of small plateaus 
(Foehn effect) formed after each step, observed in several 
magnetic molecules. Interestingly enough, we arrived at 
a convenient set of equations where the effects of dissipa- 
tion and LZS transitions can be treated separately. These 
equations account nicely for the description of both the 
magnetization steps and the plateaus. Moreover, the role 
of temperature and the field sweep rate on these effects 
becomes transparent. Finally, although our analysis is 
limited to high enough temperatures (T > 1 K) so that 
one can assume that the phonons are in equilibrium at 
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all experimental times, it nevertheless indicates how an 
enhanced Foehn effect could arise at lower T, where the 
PB effect takes place. In such a case, steep extrema in 
the M vs B curves could be observed. 

The present theoretical work provides a first step to- 
wards exploiting the possibilities that are offered by prob- 
ing magnetic molecules using external magnetic fields 
with high sweep rates. These probes, apart for provid- 
ing information specific to magnetic molecules, offer the 
possibility of conducting a detailed study of the relax- 
ational behavior of interacting spin systems as a result of 
their coupling with a "heat bath" and in particular the 
excitations of the host lattice. Development of a broad 



theoretical framework for dealing with relaxational phe- 
nomena induced by dynamical magnetic fields is indeed 
a worthy goal. 
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It is noted that the analysis of Refll5l which is based 
on a quantum master equation previously developed for 
static fields, provides an alternative set of phenomenolog- 
ical equations, to be contrasted to Eqs. I)4fi|l and Ij48ft that 
we derived from first-principles. 

This behavior should be contrasted with the sweep rate de- 
pendence of the magnetization steps for the pure quantum- 
mechanical LZS model but in the non-adiabatic regime, 
which is the relevant case for SMM's. Acording to the 
well known formula Plzs = 1 — exp[— 7r 8 2 /{2h 2 1 r)\ for the 



transition probability Plzs, when lowering r one increases 
Plzs i.e., reduces the effect of non-adiabatic transitions, 
thus giving rise to a larger height of the magnetization 
step. 

The reduction of Mout is due to thermal dissipation before 
entering the LZS regime (which determines the value of 
\Mm\ and consequently M ou t through Eq. (560)) as well 
as inside the LZS regime (shaded area in Fig. 4). 



